Nonequilibrium spin transport on Au(lll) surfaces 
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The well-known experimentally observed sp-derived Au(lll) Shockley surface states with Rashba 
spin splitting are perfectly fit by an effective tight-binding model, considering a two-dimensional 
hexagonal lattice with pz-orbital and nearest neighbor hopping only. The extracted realistic band 
parameters are then imported to perform the Landauer-Keldysh formalism to calculate nonequi- 
librium spin transport in a two-terminal setup sandwiching a Au(lll) surface channel. Obtained 
results show strong spin density on the Au(lll) surface and demonstrate (i) intrinsic spin-Hall ef- 
fect, (ii) current-induced spin polarization, and (iii) Rashba spin precession, all of which have been 
experimentally observed in semiconductor heterostructures, but not in metallic surface states. We 
therefore urge experiments in the latter for these spin phenomena. 
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I. INTRODUCTION 

Two-dimensional electron gas (2DEG) is known 
to exist in various systems, including semiconductor 
heterostructures 1 and metallic surface states. 2 Due to 
the lack of inversion symmetry introduced by the in- 
terface or surface, the spin degeneracy, the combining 
consequence of the time reversal symmetry (Kramers de- 
generacy) and the inversion symmetry, is removed and 
the energy dispersion becomes spin-split. In semiconduc- 
tor heterostructures, one of the underlying mechanisms 
leading to such spin splitting is known as the Rashba 
spin-orbit coupling, ! which stimulates a series of discus- 
sion on plenty of intriguing spin-dependent phenomena. 
Well studied phenomena include spin precession, 4,5 spin- 
Hall effect (SHE), 6 ' 7 ' 8 and current-induced spin polar- 
ization (CISP),'' 8 ' 9 ' 1() all of which have been experimen- 
tally observed in semiconductor heterostructures. Con- 
trarily, none of these in metallic surface states is reported, 
even though the Rashba effect has been shown to exist 
therein. 11 ' 12 

To the lowest order in the inplane wave vector fcy , the 
two spin-split energy branches are expressed as E± = 
h 2 k^/2m* ±afc|| (m* the electron effective mass), so that 
the Rashba spin splitting AE = E + — £L = 2akn is 
linear in fc|| . Here the proportional constant a is com- 
monly referred to as the Rashba coupling constant or 
the Rashba parameter. Typical values of a in semicon- 
ductor heterostructures are at most of the order of 10~ 2 
eV A, 7 ' 10 ' 13 while in metallic surface states a can be one 
or two orders larger. 

The first evidence of spin splitting in metallic sur- 
face states was pioneered by LaShcll et al. on 
Au(lll) surfaces at room temperature. 11 The origin 
of their observed spin splitting was later recognized 
as the Rashba effect by performing the first-principles 
electronic-structure and photoemission calculations, 1 1,1 ' 
which are in good agreement with the spin-resolved pho- 
toemission experiments. 15 ' 16 Concluded Rashba parame- 



ter of the Au(lll) surface states is about a — 0.36 eV A. 
Subsequent findings of giant Rashba spin-orbit coupling 
is also claimed in Bi(lll) surfaces 1 ' with a w 0.83 eVA 
and in Bi/Ag(lll) surface alloy 18 with a w 3.05 eV A. 

It is therefore legitimate to expect the previously men- 
tioned spin-dependent phenomena to be observed on 
those metallic surfaces with strong Rashba coupling. In 
this paper we theoretically study nonequilibrium spin 
transport in 2DEG held by Au(lll) surface states, which 
exhibit not only strong Rashba coupling but also sim- 
ple parabola- like dispersions. 19 The latter characteristic 
enables successful description of the band structure us- 
ing the simplest tight-binding model (TBM), which then 
provides the Landauer-Keldysh formalism (LKF) 20 ' 21 ' 22 
with reasonable or even realistic band parameters. 

This paper is organized as follows. In Sec. II we de- 
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FIG. 1: (Color online) (a) Tight-binding energy dispersion 
E± BM and the experimentally measured binding energy E^ p 
of Ref. 11, along the FM direction. The surface Brillouin 
zone is sketched in the inset, (b) Total density of states and 
El BM along f K and KM directions. 
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scribe the Au(lll) surface band structure by using an 
effective TBM, through which the experimentally mea- 
sured energy dispersions 11 can be perfectly reproduced. 
Section III is devoted to nonequilibrium spin transport 
on a finite Au(lll) surface channel attached to two exter- 
nal leads, using the LKF with band parameters extracted 
in Sec. II. The intrinsic SHE. the CISP, and the Rashba 
spin precession will be shown by directly imaging the lo- 
cal spin densities. We conclude in Sec. IV. 



II. Au(lll) SURFACE BAND STRUCTURE 

A. Effective tight-binding model 

We first demonstrate that the sp-derived Shockley sur- 
face states on Au(lll) from Ref. If can be well described 
by an effective TBM [sec Fig. 1(a)] for a single sheet of 
two-dimensional hexagonal lattice, taking into account 
only p z -orbital hopping between nearest neighbors, sub- 
ject to the Rashba spin-orbit coupling. The Hamiltonian 
matrix can be written as 2 5,21 



Epl 



E 1 



ikii 1 1 



V R e z ■ (ct x tj)] , (1) 



where 1 is the 2x2 identity matrix, E p is the p-orbital 
energy, tj represents the six nearest neighbor hopping 
vectors, V pp7T is the band parameter describing the orbital 
integral under the two-center approximation of Slater and 
Koster," ' Vr is the Rashba hopping parameter, and a = 
{a x , u y ,a z ) is the Pauli matrix vector. The three terms in 
Eq. (1) are the energy band offset, the kinetic hopping, 
and the Rashba hopping, respectively. Arranging the two 
primitive translation vectors for the hexagonal lattice as 
ti = (V3/2, 1/2, 0)o and t 2 = (-V3/2, 1/2, 0)o where a 
is the lattice constant, the six nearest neighbor hopping 
vectors are t/ = ±ti, ±t2, ±(ti +t2), and Eq. (1) then 
takes the explicit form of 

F*(k,|) E p + G(k ll )) 

with 

F(k||) =iVR[(l + \/3i)smk|| - ti 

+ (1 — V3i) sin k|| • t 2 + 2 sin k y a] (3) 

G(k,| ) = 2V mm [2 cos cos *j£ + cos (k y a)} . (4) 

Equation (2) can be diagonalizcd to yield the energy dis- 
persions 



£(k||)=£; p + G(k||)±|F(k||)|. 



(5) 



Noting from Eq. (3) that Vr is embedded in F(k^), the 
above dispersion contains the Rashba term to all (odd) 



In the vicinity of f, i.e., |ky|a <C 1, Eqs. (3) and 
(4) are approximated by F(k||) ss — 3Vr (k x — ik y ) a and 
G(k||) f=a 6V pp7r — {ZVppnCL 2 /2)kp respectively, and the 
Hamiltonian matrix (2) then takes the form 



E 



3V, 



-3V R {k x - iky) a 



-W R (k x +ik y )a E a 



3Vd rj7r^ 



(G) 



where Eq = E p + 6V pp7T . Equation (6) is consistent 



with the p z -resolved effective Hamiltonian of the earlier 
TBM by Petersen and Hedegard, who considered all the 
three p-orbitals, subject to the intra-atomic spin-orbit 
coupling. 21. 



B. Extraction of band parameters 

We now fit our tight-binding dispersions (5) with the 
experiment of Ref. ] ! . This can be done by comparing 
the low-fc|| expansion of Eq. (5), 



£( k ll)| 



fea<Cl 



E p + 6V m 



3 Vrm-TrQ 



fcf±3Vi?a|k|||, (7) 



orders in k\\ . 



with that of the free-electron model, E(k\\) = Eo + 
(h 2 /2m*)k^ ± ak\\. In addition to the band offset Eq 

/'.',, + GVppx, we identity V pp7r = — (2/3a 2 )(S 2 /2m*) 
and Vr = a /3a. Using the reciprocal vector gi = 

(4 7 r/v / 3a)(l/2,v / 3/2,0) and M = 1.26 A -1 from Ref. 

11, we have h 2 /2m* w 15.2 eVA 2 , a « 0.3557_eVA, 
and Eq ~ —0.415 eV. The norm of gi gives M such 
that the lattice constant is a = 47r/\/3|gi| = 5.7581 A. 
Hence we deduce V pp7T = -0.3056 eV, Vr = 0.0206 eV, 
and E p = 1.4188 eV. Substituting these parameters into 
Eqs. (3)-(5), a nearly perfect consistency between our 
effective TBM and the experimentally measured binding 
energy of Ref. 11 can be seen in Fig. 1(a). The ex- 
perimentally measured Fermi surface of the concentric 
rings slightly distorted from circles 19 can be reproduced 
as well, but we do not explicitly show. 



III. NONEQUILIBRIUM SPIN TRANSPORT 

A. Landauer-Keldysh formalism vs tight-binding 
model 

Next we apply the Landauer-Keldysh formalism, 21 
namely the nonequilibrium Keldysh Green's function 
formalism - ' applied on Landauer multitcrminal ballistic 
nanostructures. For detailed introduction to the LKF, 
see Refs. 20,22. To make use of the previously extracted 
band parameters in the LKF calculation, we consider the 
second-quantized single particle Hamiltonian, 21 

U = ^2 Sncic,, + ^2 c) m [t + Hr (ct x d lnn ) J c„, (8) 

n (m,n) 
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FIG. 2: (Color online) Local spin density of (a) the out- 
of-plane component (S z ) and (b) the inplane component 
((5a,), (S y )) in a 398.9 X63.3 (total number of sites N = 931) 
conducting sample made of Au(lll) surface. The size of 
each local marker depicts the magnitude. In (a), red/dark 
(green/light) dots denote (S z ) > ((&) < 0). The maxi- 
mum of (S z ) is 1.2 x 1CT 3 (ft/2) while the mean of (S y ) is 
3.16 x 10~ 4 (ft/2). 



which is equivalent to Eq. (1), provided e n = E p , 
to = ~Vpp n , and t R = V R . In Eq. (8), cj, (c„) is the 
creation (annihilation) operator of the electron on site n, 
(m, n) means that sites to and n are nearest neighbors to 
each other, and d mn is the unit vector pointing from n to 
to. Despite the different system sizes TBM and LKF con- 
sider (infinite for TBM and finite for LKF) and different 
functions they provide (simple band calculation by TBM 
and nonequilibrium transport by LKF), the equivalence 
of the underlying Hamiltonians should contain the same 
physics. The explicit correspondence can be shown by 
comparing the band structure by the TBM with the to- 
tal density of states (TDOS) by the LKF, provided that 
the same parameters are used. 

For the LKF calculation, we consider a 80(a\/3/2) x 
11a « 398.9 A x63.3 A (total number of sites N = 931) 
channel made of an ideal Au(lll) surface, in perfect con- 
tact with two unbiased normal metal leads at the left 
and right ends of the sample. We will further image the 
nonequilibrium spin transport on this two-terminal setup 
later. As shown in Fig. 1(b), the range of the calculated 
nonvanishing TDOS is consistent with the TBM disper- 
sion along the TK direction, which corresponds to the 
nearest neighbor hopping direction as we considered in 
the underlying Hamiltonian (8). 



B. Injection of unpolarized current: Intrinsic 
spin-Hall effect and current-induced spin 
polarization 

Combination of the consistency between the experi- 
mental and the TBM dispersions, and that between the 
dispersion by the TBM and the TDOS by the LKF, in- 
directly demonstrates that the following imaging of lo- 
cal spin densities by the LKF stands on an experimental 
footing. As a first demonstration of the nonequilibrium 
spin transport, we turn on the bias of potential differ- 
ence eVb = 0.2 eV between the two normal metal leads. 
We will denote ±eVb/2 on the leads by ± sign. With 



such injection of an unpolarized electron current, we ex- 
pect (i) the SHE of the intrinsic type, and (ii) the CISP, 
which follows the Rashba eigenspin direction of the lower 
energy branch. Both of these can be seen respectively 
in Figs. 2(a) and 2(b). The former shows an antisym- 
metric out-of-plane spin accumulation with (S z } max = 
1.2 x 10~ 3 (fi./2) at lateral edges, while the latter shows 
that the inplane components of spins mostly point to 
+y axis 28 with average value (S y ) = 3.16 x 10~ 4 (fr/2). 
Note that the local spin densities shown here represent, 
by definition, the site-dependent total number of spins. 
Dividing (S y ) by the hexagonal unit cell area v / 3a 2 /2, 
we deduce that the obtained spin (area) density due to 
CISP is in average 1101 /im~ 2 , which is clearly much 
stronger than that observed in the CISP experiment of 
Rcf. 9, where the spin (volume) density less than 10 
/im~ 3 (corresponding to a even weaker spin area density) 
is reported. 



C. Injection of spin-polarized current: Rashba spin 
precession 

Next we inject spin-polarized currents by replacing the 
left (source) lead with a ferromagnetic electrode. Previ- 
ously, the self-energy due to the normal metal lead, which 
is assumed to be semi-infinite, in thermal equilibrium, 
and in perfect contact with the sample, can be obtained 
by solving the surface Green's function of the lead, 20 sub- 
ject to Hamiltonian Hi c &d = P 2 /2m + V. The momentum 
operator p = (p x ,p y ) is two-dimensional and the poten- 
tial V describes a infinite potential well of a semi-infinite 
rectangle shape. The exact form of the lead self-energy 
reads 



£ fl (pi, P2 ) 



2t 



nirpi . nirp2 
sin sin ■ 



N d + 1 ^ N d + 1 N d + 1 



X 1 



glfc„a gm (fc raC j) 

k n o, 



(9) 



with k n a = yj (E - E n ) /t and E n = [mr/(N d + l)]H. 
Here Pi(2) = E 2, • • • , N d is the lateral position (implic- 
itly in units of lattice constant a) of the edge site in the 
sample in contact with the lead, t is the coupling be- 
tween the sample and the lead and is usually set equal 
to the kinetic hopping to in the sample, and (Nd + 1) a 
is implicitly assumed to be the width of the lead. 

To take into account the exchange field inside the fer- 
romagnetic lead, we adopt the Weiss mean field approx- 
imation and add a Zeeman term —[Xb<? • B ea; (hb ~ 
5.8 x 10- 5 eVT _1 the Bohr magneton) to Hi oa d- Typ- 
ical exchange field may be as high as |B ea; | ps 10 -3 T 
(Ref. 29), leading to |^ B B e:c | ps 5.8 x 10~ 2 cV. Wc 
will take this value in the forthcoming spin precession 
demonstration. The explicit form of the self-energy is 
obtained by substituting in Eq. (9) k n a — -> k°a = 
y/(E-E n + a\ti B B ex \)/t and 1 -► £ ff=± \mr) (na\, 
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FIG. 3: (Color online) Local spin density of (a) the out- 
of-plane component (S z ) and (b) the inplane component 
((St)) (Sy)), in a conducting sample made of Au(lll) surface, 
subject to a ferromagnetic source lead with +x magnetization, 
(c) (S x ), {Sy), and (S z ) as a function of x at y = 0, i.e., along 
the dashed line sketched in (a). Computed values are com- 
pared with the previously obtained spin vector formula based 
on quantum mechanics. 



where |ner) is the spin- 1/2 state ket with quantization 
axis n. 30 

Applying the same voltage difference of 0.2 V and mag- 
netizing the ferromagnetic source lead along +x axis, 
Figs. 3(a) and 3(b) show the out-of-plane and inplane 
components of the local spin densities, respectively. The 
injected x-polarized spins moving along +x and encoun- 
tering the Rashba effective magnetic field pointing to —y, 
are forced to precess about — y-axis counterclockwise, and 
hence the Rashba spin precession is observed. In the free 
electron model, the spin precession length (the distance 
within which the spin completes a precession angle of tt) 
is L so = (ir/a)(h 2 /2m*) w 134 A. Thus the channel 
length 398.9 A is about 3 times L so , which is consistent 
to what we observe in Figs. 3(a) and 3(b). Note that here 
the SHE competing with the spin precession is relatively 
weak due to the strong exchange field we consider in the 
source lead. However one can still observe the tiny asym- 
metry of the (S z ) pattern of Fig. 3(a) along the lateral 
direction (more +(S Z ) and —(S z ) accumulations near the 
bottom and top edges, respectively). In the following we 
will concentrate on the spin precession only. 

To compare the LKF results with the free electron 
model in further detail, we recall the spin vector for- 
mula [see Eq. (6) of Ref. 31], which takes the form of 
({S x ),{S y ),(S z )) = (cosA0,O,sinA0) here with A9 = 
(2m*/h 2 )ax = 2.34 x 10~ 2 x, where x is in unit of A. 



Note that a factor of \/3/2 responsible for the net and 
actual hopping distances has to be taken into account in 
x, since the crystal structure information remains. Ac- 
cordingly, good agreement between the LKF and the spin 
vector formula can be seen in Fig. 3(c). Note that in 
view of both Figs. 2 and 3, size and edge effects, arising 
from the charge distribution, are also observed. The for- 
mer, the size effect, appears in the modulation along y 
direction with roughly 4 peaks corresponding to a wave 
lengths 32 A, roughly shorter than the Fermi wavelength 
27r/fcF « 38 A; the latter, the edge effect, appears in the 
abnormal charge accumulation near the side and drain 
edges. 



IV. CONCLUSION 

In conclusion, we have shown that the sp-derived 
Shockley surface states on Au(lll), 11 ' 14 ' 15 ' 16 ' 19 which ex- 
tend over the first few layers though, can be well de- 
scribed by an effective TBM for a two-dimensional hexag- 
onal lattice, taking into account p z -orbital and near- 
est neighbor hopping only. Required parameters in the 
nonequilibrium spin transport calculation by the LKF, 
demonstrating (i) intrinsic SHE and CISP due to injec- 
tion of unpolarizcd current and (ii) the Rashba spin pre- 
cession due to injection of spin-polarized current, thus 
stand on an experimental footing of the pioneering work 
of LaShell et. al.. Ll Calculated local spin densities in 
all the three spin phenomena arc much stronger than 
those in semiconductor heterostructures. Whereas the 
magnetic optical Kerr effect (MOKE) can sensitively de- 
tect a spin volume density of less than 10 spins per /im 3 
(Ref. 9), our results of more than 10 3 spins per /im 2 sug- 
gest definitely measurable nonequilibrium spin transport 
supported by the Au(lll) surface states and others with 
even stronger Rashba coupling such as Bi(lll) surfaces 1, 
or Bi/Ag(lll) surface alloys. lf " Last, in addition to the 
stronger local spin densities induced by stronger Rashba 
coupling, the spin precession length (typically of the or- 
der of 1 /im in semiconductor heterostructures) in these 
surface states is greatly reduced [134 A for Au(lll) re- 
ported here] , such that the fine structure of the spin pat- 
terns due to spin precession or the intrinsic SHE requires 
high resolution apparatus such as the spin-polarized scan- 
ning tunneling microscopy. '" 
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